Fractal Dynamics in Chaotic Quantum Transport 
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Despite several experiments on chaotic quantum transport, corresponding ab initio quantum sim- 
ulations have been out of reach so far. Here we carry out quantum transport calculations in real 
space and real time for a two-dimensional stadium cavity that shows chaotic dynamics. Applying 
a large set of magnetic fields yields a complete picture of the magnetoconductance that indicates 
fractal scaling on intermediate time scales. Two methods that originate from different fields of 
physics are used to analyze the scaling exponent and the fractal dimension. They lead to consistent 
results that, in turn, qualitatively agree with the previous experimental data. 
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Since the pioneering works of Mandelbrot [1], fractal 
patterns have been found in a variety of objects in nature 
including, e.g., snowflakes, fern leaves, coastlines [2j|3], 
and even music [4]-[7]. These self-similar (or self-affine) 
structures were also found in many branches of chem- 
istry and physics; prominent examples are crystal growth 
and fractal surfaces, and transport in gold nanowires and 
electron "billiards" [3j I8HTT]. In contrast with ideal- 
ized mathematical fractals continuing to infinitely small 
scales, fractal scaling in nature has a lower and an upper 
limit. 

While fractals found in nature are often well described 
by classical theories [U [3] [BJ [9] , fractals have also been 
suggested to manifest in different quantum systems [T2l - 
[17], where a fundamental lower cutoff for fractal scaling 
is given by the Heisenberg uncertainty principle. In the 
case of transport through chaotic systems, such as chaotic 
electron billiards, both semiclassical (involving quantum 
interference) and classical mechanisms for the emergence 
of fractal scaling have been proposed p~8j [19] . 

For quantum systems with an underlying classically 
mixed phase space with both regular and chaotic regions, 
a quantum graph model suggests a splitting of the chaotic 
regime into two parts [14 : one part yields fractal conduc- 
tance fluctuations while the other one leads to isolated 
resonances on small scales. These isolated resonances 
were later shown to be associated with the eigenstates of 
a closed system [20] . 

A stadium quantum billiard of charged particles is a 
generic chaotic system, whose underlying classical phase 
space is chaotic. The phase space becomes mixed in pres- 
ence of a (perpendicular) magnetic field. In the past 
two decades the system has been subject to several ex- 
periments [10] [21-24 . A typical setup consists of the 
two-dimensional electron gas (2DEG) in a semiconduc- 
tor heterostructure, where metallic gates are used to form 
the geometrical shape of the "billiards" - here called a 
quantum dot. Alternatively, stadium billiards (and other 
chaotic systems) can be realized experimentally with mi- 



crowave cavities [24] . 

In this Letter we calculate the fractal scaling of conduc- 
tance fluctuations in an open quantum stadium billiard. 
Our explicit solution of the time-dependent Schrodinger 
equation in real space and real time for chaotic transport 
goes beyond both the semiclassical treatment [18 and 
the above mentioned quantum graph model [14] . We an- 
alyze the fractal scaling using a method which we call 
minimum- maximum (MIN-MAX) method [3j [10] [25] , as 
well as detrended fluctuation analysis [26-28J (DFA). The 
MIN-MAX method was used by Sachrajda et al. [10] for 
the analysis of experimental magnetoconductance curves. 
We are able to find a good agreement between the- 
ory and experiment, both yielding a fractal dimension 
D~ 1.2... 1.3. 

Model and the computational scheme.- We consider 
a model for semiconductor stadium device fabricated in 
the 2DEG of a AlGaAs/GaAs heterostructure similar to 
Ref. [10 . The Hamiltonian describing our 2D system 
reads (in atomic units) 

H=\[-iV + A(r)} 2 + V ext (r,t), (1) 

where the vector potential is given in the linear gauge, 
A(r) = {—By, 0,0), to describe a static and uniform 
magnetic field perpendicular to the plane. During the 
time-propagation at t > 0, the potential V ex t(^^t) con- 
sists of three parts: (i) a stadium with radius r = 1 
and width d = 0.7, (ii) input and output leads of width 
w = 0.56, and (iii) a linear potential along the propa- 
gation direction describing a source-drain voltage. The 
potential has hard boundaries with a depth Vb = 10000 
and the slope of the accelerating linear potential is —100. 
The central part of the external potential is shown in 
Fig. [I] The input and output lead extend further to the 
left and right. 

The initial state at t = is calculated by taking a small 
part of the input lead as a potential well. The resulting 
ground state of a single electron in the well is then used 
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FIG. 1: Snapshots of the electron density in the model sta- 
dium system (see text) during a transport simulation with the 
magnetic flux $/$o = 20. The input and output leads extend 
further to the left and right. 



as an initial state for the time propagation. At t > 
the above described linear potential accelerates the wave 
packet across the system. For the time propagation we 
use a fourth-order Taylor expansion of the time-evolution 
operator. The octopus code package [29 - previously 
applied for similar real-space transport calculations in 
quantum rings [30] and Aharonov-Bohm oscillations [3TJ 
[32] - is used in all the calculations. 

We assess the conductance by calculating the inte- 
grated probability density in the output lead from 



T($,t) 



di#($,M)| 5 



(2) 



output 



where is the fixed magnetic flux, given above in units 
of the magnetic flux quantum <l>o = h/e. We call T as 
a transmission coefficient, whose validity in estimating 
the relative conductivity as a function of an external pa- 
rameter - here the magnetic flux - has been justified in 
Ref. [30 . Thus, we repeat the time-propagation for dif- 
ferent values of <E> to obtain the magnetoconductance that 
can be compared with the experiments in Ref. [TO] , 

Transport simulations.- In Fig. [T] we show snapshots 
of the electron density at different times at $ = 20$o 
through the stadium. Approximately one half of the 
density is transferred through and other half is either re- 
flected back to the input lead or confined in the stadium. 
As expected, the density is scattered in the stadium in 
a chaotic fashion. The size of the wiggles during the 
scattering depends on the momentum - the higher the 
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FIG. 2: Transmission coefficient as a function of time and 
magnetic flux through the stadium. 



momentum the higher eigenstates are probed. Here, we 
have chosen the momentum such that sufficient qualita- 
tive complexity is observed, which, as shown below, leads 
also to a complex behavior of T. On the other hand, the 
momentum is limited by the grid spacing of the simu- 
lation box - all the nodes in the scattered wave packet 
should be accurately described. 

A complete presentation of our transport results is 
given in Fig. [2] where the transmission coefficient T is 
plotted as a function of both time and the magnetic flux. 
The figure consists of 401 respective time-propagations, 
each with a fixed number of flux quanta ^/^o ranging 
from zero to 40 in steps of 0.1. The flux range is qualita- 
tively similar to the experiment in Ref. [10]. A complex 
magnetoconductance is formed if the propagation time 
is larger than ~ 1. A cross section of the conductance 
at t = 1.4 is shown in Fig. [3] We point out that due to 
the finite system size we are not able to reach the equi- 
librium current and thus find the absolute conductance. 
Therefore, we consider fixed propagation times through 
the parameter range of ^/^o- In other words, a fixed 
propagation time is expected to treat all the values of 
^/^o equally in order to obtain the relative conductance 
T. 

Now, the essential question is after which times fractals 
are found, and how the fractal dimension changes over 
time. This is analyzed in the following with two tech- 
niques: the MIN-MAX method [3] [TO] [25], and DFA [26]- 

Minimum-maximum method.- To extract the fractal 
dimension D for a mapping / : M R, the domain 
of the given function is first divided into length intervals 
Ax. The difference between the minimum and maximum 
of the function is calculated within every interval and 
added up. In case of fractal scaling, the resulting sum is 
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FIG. 3: Transmission coefficient as a function of the magnetic 
flux at t — 1.4. The inset shows the scaling exponent a — 1.46 
obtained from the DFA analysis. 



a power-law function of the interval length [3j [TO] [25 
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Detrended fluctuation analysis.- DFA is a standard 
method that was developed in the context of time-series 
analysis [26 and has proven to be very reliable par- 
ticularly in dealing with non-stationary time-series and 
trends in the data [6] [26] EH [33j [34] . It has also been used 
outside the time domain, e.g., to study the organization 
of DNA nucleotides [28]. However, to our knowledge, 
DFA has not been applied to fractal conductance curves 
before. 

The standard procedure of DFA consists of the follow- 
ing four steps [26] [27]: (1) integrating the time series, (2) 
dividing the series into windows of size s, (3) fitting with 
a polynomial f s (i) of degree m = 2 . . . 4 that represents 
the trend in the window, and (4) calculating the variance 
with respect to the local trend f s (i) from 



F(s) = ((f(i)-fs(i)) 2 ) 



1 N 

2=1 



(4) 



The key point in applying DFA to study conductance 
fluctuations is to relate the exponent a to the quantity 
of interest (here: the fractal dimension D). It is known 
that D = 2 - 7 /2 with ((AG) 2 ) - (AB)i [TO] [18]. The 
latter is exactly step (4) of the DFA analysis above. We 
therefore omit step (1) and identify a = 7, hence the 
fractal dimension reads D = 2 — a/ 2. 

Results.- In DFA, we apply quadratic detrending 
(m = 2) to our data in Fig. [3J The inset shows the 
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FIG. 4: Fractal dimension D calculated from DFA with the 
relation D = 2 - a/2 and from the MIN-MAX method, re- 
spectively, during the time-propagation. Note that the frac- 
tal structure is developed only at t > 1. The inset shows 
the time-development of the error in the fitting procedure at 
t = 0.8. . . 1.5 (see text). 



fitting of the data (solid line) at t = 1.4 that yields 
a = 1.46. This qualitatively agrees well with the experi- 
mental result 7 = a = 1.38 of Sachrajda et al. [10 . The 
corresponding fractal dimension extracted from DFA is 
D = 1.27. In comparison, the MIN-MAX method yields 
D = 1.32 for our data, whereas the corresponding ex- 
perimental result - obtained with the same method - is 
D = 1.25 [TO]. Concluding, we find an excellent agree- 
ment of the results both regarding the different methods 
and comparison with the experimental data. We point 
out that our stadium model is similar to the experiment 
and the channel dimensions are also comparable. Ac- 
cording to our calculations, increasing the channel width 
from 0.56 to 0.7 leads to the same D obtained in the MIN- 
MAX method, whereas DFA yields a slightly smaller D. 

In Fig. [4] we show the time-development of the frac- 
tal dimension obtained from DFA and the MIN-MAX 
method, respectively. We point out that at times t < 0.5 
the quality of the fitting is poor - or the fitting is ill- 
defined per se - as the conductance curve as a function 
of the flux has not yet been developed (see Fig. |2|. Sec- 
ondly, the clear signatures of a fractal structure are de- 
veloped only at t > 1. Nevertheless, D converges dur- 
ing the time-propagation towards the values given above, 
and the quality of the fitting in both methods improves 
as well. The inset in Fig. [4] shows the error in the fit- 
ting, ERR = 1 - CC 2 , for times t = 0.8 . . . 1.5. Here CC 
is the Pearson product-moment correlation coefficient of 
the log-log data. Thus, ERR measures the linear fit qual- 
ity such that ERR = corresponds to exact linear be- 
havior. The minimum of the error is obtained at t « 1.4, 
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which is the optimal time used above to determine a and 
D. At larger times with 1.45 < t < 1.5, we find a small 
bump in D in Fig. [4] corresponding to an upward kink 
in the error in the inset. This is the first signature of 
back-scattering effects due to the finite simulation box 
(see above). In this way we are able to determine the 
range of validity in our scheme to determine the fractal 
dimension. 

In summary, we have calculated the time-evolution of 
a single-electron wave packet through a two-dimensional 
stadium-shaped cavity by solving the Schrodinger equa- 
tion in real time and real space. The relative conductance 
has been calculated for a large set of magnetic fluxes in 
order to analyze the fractal nature of the magnetoconduc- 
tance. We have found that the conductance shows clear 
indications for fractal scaling. The fractal dimensions ex- 
tracted from two respective methods are consistent with 
each other. Moreover, we have found an excellent quali- 
tative agreement with previous experimental results. Our 
findings indicate that DFA suits well for the analysis of 
fractal scaling in chaotic quantum transport. Hence, we 
suggest to extend the use of the concept of data defend- 
ing (and hence DFA) to study fractal scaling of transport 
and other characteristics in chaotic (quantum) systems. 
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